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ABSTRACT 

Rigorous  expressions  for  diffusion-controlled  currents  at  a  sta¬ 
tionary  finite  disk  electrode  are  derived  through  use  of  the  Winer-Hofp 
technique.  The  chronoamperometric  curve  obtained  varies  smoothly  from  a 
curve  represented  by  the  Cottrell  equation  to  a  steady  state  value 
similar  to  that  obtained  for  a  spherical  electrode  as  time  elapses.  The 
solution  can  be  expressed  also  as  the  Cottrell  term  multiplied  by  a 
power  series  in  the  parameter  *T)t/a,  where  a  is  the  electrode  radius. 

The  present  work  is  discussed  in  terms  of  the  coefficient  of  the  first 
term  in  this  formulation. 
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INTRODUCT ION 

Most  studies  of  mass  transfer  in  electrochemical  problems  have 
focused  on  linear  diffusion  to  and  from  a  planar  electrode  or  on  diffusion 
in  the  radial  direction  at  a  spherical  or  a  cylindrical  electrode.  In 
other  words,  diffusion  is  considered  only  in  the  direction  perpendicular 
to  an  electrode,  though  the  electrodes  employed  in  practice  have  their 
own  more  complex  geometries.  Real  electrode  processes  are  frequently 
complicated  by  diffusion  with  directional  components  other  than  normal 
to  the  electrode,  i.e.,  non-linear  diffusion,  especially  near  the  edges 
of  the  electrode.  We  can  take  as  examples  of  complicated  diffusion 
shielding  effects  at  the  DME  or  the  HMDE  by  the  glass  tip  supporting  the 
mercury  drop  [1-3],  edge  effects  at  planar  electrodes  [4-15]  and  effects 
of  partially  blocked  electrodes  [16-24].  These  effects  are  strongly 
dependent  not  only  on  geometries  of  electrodes  but  also  on  the  relation 
between  dimensions  of  electrodes  and  the  time  scale  of  the  experiment. 

One  of  the  simplest  geometries  of  electrodes  that  complicates 
diffusion  is  a  planar  disk  electrode  embedded  flush  in  an  infinite  in¬ 
sulating  plane.  Chronoamperometric  curves  at  this  electrode  are  ex¬ 
pected  to  have  the  following  behavior  as  time  elapses:  at  sufficiently 
short  times  they  should  obey  the  Cottrell  equation  (i  'vl/v’F-),  be 
gradually  influenced  by  three  dimensional  diffusion  near  edges  of  the 
electrode  at  longer  times,  and  finally  approach  a  steady  state  current, 
as  is  observed  at  a  spherical  electrode  [25].  One  of  the  interesting 
points  of  the  stationary  disk  electrode  Is  that  the  steady  state  current 
can  be  observed  at  long  times  though  the  electrode  is  stationary  as  well 
as  planar.  Another  interesting  point  is  that  planar  disk  electrodes  can 
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be  constructed  easily  to  conform  closely  to  a  given  model.  This  is  in 
contrast  to  hanging  mercury  drop  electrodes  or  shielded  planar  elec¬ 
trodes.  The  other  fascinating  point  is  that  for  proper  dimensions  and 
times  such  electrodes  make  it  possible  to  determine  n  and  D  simultaneous¬ 
ly  from  a  single  chronoamperogram  [9,10,28,29]. 

Both  theoretical  and  experimental  work  has  been  carried  out  on  mass 
transport  at  the  stationary  disk  electrode.  Bard  [4]  concluded  from 
potentiometric  measurements  that  a  planar  disk  electrode  has  a  transi¬ 
tion  time  which  is  similar  to  that  obtained  at  a  spherical  electrode 
when  the  time  scale  is  long  enough.  Lingane  [5]  gave  a  quantitative 
experimental  description  of  the  chronopotentiometric  and  chronoampero- 
metric  constants  at  unshielded  planar  disk  electrodes.  Soos  and  Lingane 
[6]  derived  an  analytical  expression  for  diffusion-1 imited  currents  at  a 
planar  disk  electrode.  However,  their  equation  is  just  the  sum  of  the 
Cottrell  equation  and  the  expression  for  the  steady  state  current,  so 
that  it  does  not  describe  actual  chronoamperometric  curves.  Flanagan 
and  Marcoux  [7]  evaluated  departure  of  chronoamperometric  curves  from 
the  Cottrell  term  by  means  of  digital  simulation.  Ito  et  al.  £8]  in¬ 
troduced  an  empirical  parameter  into  an  expression  for  the  diffusion 
current  as  a  result  of  chronoamperometric  experiments  at  platinum 
microelectrodes.  Kakihana  et  al.  [9,10]  did  potentiostatic  experiments 
at  unshielded  small  disk  electrodes  in  order  to  obtain  more  accurate 
values  of  diffusion  coefficients  of  electroactive  species.  They  also 
employed  digital  simulation  using  a  corrected  version  of  the  scheme  of 
Flanagan  and  Marcoux  [7]  and  compared  their  simulated  values  with  their 
experimental  ones.  Recently  Sarangapani  and  DeLevie  [11]  challenged 
theoretically  the  problem  of  ac  response  at  this  kind  of  electrode; 
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their  expression  does  not  satisfy  the  boundary  condition  given.  Dayton 
et  al.  [12]  measured  almost  time  independent  diffusion  current  at  small 
sizes  of  carbon  fiber  electrodes. 

Most  recently  Oldham  has  analyzed  theoretically  the  problem  of  current 
distribution  for  an  infinitesimally  thick  electrode  bounded  in  only  one  of 
the  remaining  four  directions  [14].  The  solution  can  be  applied  through 
appropriate  transformation  to  finite  disk  electrodes  of  moderate  size.  In 
addition,  Heinze  has  developed  a  digital  simulation  model  for  current  at 
a  planar  disk  electrode  [15]. 

Theoretical  studies  of  complicated  diffusion  have  been  done  on 
assumptions  that  potential  theory  valid  only  for  the  steady  state  may  be 
applicable  to  time  dependent  systems  [17-21],  and  non-linear  parts  of 
diffusion  may  be  represented  as  mean  concentrations  [22-24].  Neither 
technique  is  a  direct  approach  to  non-linear  diffusion.  In  the  field  of 
heat  transfer,  which  is  similar  to  the  electrochemical  problem  from  the 
phenomenological  point  of  view,  a  rigorous  solution  has  been  given  to 
systems  with  intricate  geometrical  boundaries  [26]  by  employing  the 
Wiener-Hopf  method  [27].  This  technique  is  not  only  powerful  but  also 
applicable  to  many  electrochemical  problems  concerning  complicated 
diffusion  with  various  kinds  of  geometries  of  electrodes. 

The  purpose  of  this  paper  is  to  derive  rigorous  expressions  for 
diffusion-controlled  currents  at  the  unshielded  planar  disk  electrode  in 
quiescent  solution. 

DERIVATION 

Let  us  consider  a  simple,  reversible  electrode  reaction,  0  +  ne  + 

R,  involving  only  species  soluble  in  the  solution.  We  assume  that 

diffusion  coefficients  of  both  species  have  the  common  value,  D,  and  ^ 

J 


that  migration  can  be  neglected.  Under  these  conditions,  the  expres¬ 
sions  for  diffusion-controlled  currents  at  the  planar  disk  electrode  are 
derived  for  large  values  of  t  and  for  small  values  of  t. 

Large  Values  of  t  (Descending  Series).  The  diffusion  equation 
represented  in  polar  coordinates  is  of  the  form 

* 0  <7  If  <r2  if1'  4  -7—^  Is (si"6  ss1'1  (j  * 0  °r  R)  <’> 

where  C.  =  C.  (r,  9,  t)  is  the  concentration  of  species  j,  r  is  the  dis- 

J  J 

tance  from  the  center  of  the  electrode  and  6  is  the  angle  which  the 
radius  vector  makes  with  the  z-axis  as  shown  in  Fig.  1. 

The  initial  conditions  are  given  by 

CQ(r,  6,  0)  =  C°,  CR(r,  6,  0)  -  0  (2) 

where  C°  is  the  bulk  concentration  of  species  0.  Taking  into  account 
that  the  Nernst  equation  holds  at  the  electrode,  we  have  the  boundary 
condition  for  e  =  tt/2 

CQ(r,  tt/2,  t)  =  cCR(r,  tt/2,  t)  for  0  <  r  <  a  (3) 

where  z,  =  exp[(nF/RT)(E  -  EQ)]  and  a  is  the  radius  of  the  electrode. 
Since  the  sum  of  the  fluxes  for  the  two  species  at  the  electrode  is 
equal  to  zero,  it  follows  that 

3C0(r,  n/2,  t )/36  +  3CR(r,  tt/2,  t)/39  =  0  for  0  <  r  <  a  (4) 

The  boundary  conditions  at  the  insulating  plane  are  given  by 

(r,  tt/2,  t)/39  *  0  for  r  >a  (j  *  0  or  R)  (5) 

If  Jr~  (CQ  -  C°  ;  CR)  is  replaced  by  v  and  the  Laplace  transforma¬ 
tion  is  carried  out  with  respect  to  t,  eqns.  (1),  (2),  (3)  and  (5)  are 
reduced  to 


(r2s/D  ♦  1/4)  v  ■  r3  (r3v/3r)/3r  ♦  (l/sin0)3(sin63v/39)/3e 


(6) 
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v(r,  tt/2)  =  j-  /r  for  0  <  r  <  a 

(7) 

3v(r,  it/2)/30  =  0  for  r  >  a 

where  v  means  the  Laplace  transform  of  v.  If  we  apply  the  Kontorovich- 
Lebedev  transformation  [30]  given  by 

oc 

=  /  v(r ,6)K  (r/s7D)(l/r)dr  (8) 

0  y 

to  eqn.  (6),  taking  into  account  v  0  (vHF)  as  r  ■*  0,  we  have 
(y2-l/4)$  +  (1/sinQ)  (3  (sine3cp/36)/a8)  =  0  (9) 

Although  eqn.  (9)  holds  only  for  |Re  y|  <  1/2,  analytical  continuation 
[31]  makes  it  valid  for  any  values  of  y.  The  solution  of  eqn.  (9)  is 

~  B-j  +  B2  6 ) 

where  P  and  Q  are  the  first  and  the  second  kind  Legendre  functions  [32], 
respectively,  and  and  B2  are  any  constants  independent  of  6.  Since 
<{>  should  be  continuous  at  0  =  0,  B2  is  found  to  be  zero.  Eliminating  B-j 
in  terms  of  the  combination  of  $  and  3(f>/30  and  letting  0  be  tt/2  yields 

Py-l/2(°)^(u,Tr/2)  +  pu.l/2{0)^(y,1T/2)  =  0  (10) 

where  the  prime  means  the  derivative  with  respect  to  0. 

By  considering  eqns.  (7)  and  (8),  an  alternative  expression  for 


eqn.  (10)  is  written  as 


-7- 


Q  00 

-(C°/s)  I  K  (rV  s/D)  (l/vr7)dr  +  /  v(r,-n/2}  K  (r/  s/D)(l/r)dr 
0  y  a  y 

=  L(u)  f  v'  (r, tt/2)  K  (r/  s/D)(l/r)dr  (11) 

0  u 

with 

Hx)  ■  -p1.-V2(0)/PVl/2(0)  = 

r(l/4  +  u/2)r(l/4  -  u/2)/(2r(3/4  +  u/2)r(3/4  -  u/2))  M2) 

Now,  a  new  function  $  is  defined  as 

*(u)  =  rd  +v)r~u  f  v'(r,7T/2)  I  (r/17P)  (l/r)dr  (13) 

0  y 

where  E,  =  (a/2)  A/D  and  I(j  is  the  modified  Bessel  function  of  the 
fjrst  kind.  t(u)  is  found  to  be  regular  in  the  region  Re  y  >  -1/2, 
because  I  (z)  becomes  (z/2)u/r(l+u)  [32]  as  u  approaches  infinity,  and 

M 

v'  has  the  order  of  /r  as  r  approaches  zero.  On  making  use  of  eqn.  (A- 
2)  in  Appendix  1,  the  integral  on  the  right  hand  side  of  eqn.  (11)  can 
be  expressed  in  terms  of  4>  as 

/  v'  (r,  J)  Kp  (/§  r)  ^  =  (l/2)C‘wr(p)«(-u)  +  (l/2)CWr(-u)  *(w)  (14) 

l(u)  defined  by  eqn.  (12)  can  be  decomposed  into  the  product  of  two 
functions: 


L(u)  «  Lp(u)/LN(u) 


(15) 
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where 

Lp(y)  =  r(l/4  +  y/2)/2r(3/4  +  y/2) 

Ln(u)  =  r( 3/4  -  u/2)/r(l/4  -  u/2) 

Lp(u)  and  L^(p)  are  analytic  and  non-zero  in  the  regions  Re  y  >  -1/2  and 
Re  y  <  1/2,  respectively.  Substituting  eqns.  (14)  and  (15)  into  eqn. 
(11)  and  multiplying  both  sides  by  Ljj(y )/Cyr(-y )  yields 

V  +  X  =  W  +  (1/2)  Lp(y)  v(y )  (16) 

V (y )  =  ( Ln ( y )/ C^r(-y) )  /  (-C°/s)K  ,  (r/17U)(l/^r)  dr  (17) 
n  0  y 

X(y)  *  (LN(y)/^r(-y))  f  v  (r,-n/2)  K  (r/  s/D)(l/r)dr  (18) 

a  M 

W(y)  *  Lp(y)r(y)*(-y)/2<-2lir(-y)  (19) 

Since  two  unknowns,  v(r,ir/2)  and  v'(r,ir/2)  are  involved  in  eqn. 
(16),  explicit  solutions  cannot  be  obtained  in  an  ordinary  algebraic 
manner.  However,  the  Wiener-Hopf  technique  [27]  sometimes  enables  us  to 
divide  one  equation  into  two  independent  ones  leading  to  solutions.  If 
we  apply  Cauchy's  integral  theorem  [31]  to  V (y )  for  -1/2  <  Re  y  <  0  and 
take  the  imaginary  part  of  the  path  of  integration  to  infinity,  we  have 

V(y)  -  Vp(y)  +  VN(y)  (20) 

f.+i"  , 

Vp(y)  -  -(l/2iri )  /  V(z)(z-y)  'dz 

B-i°° 

VN(y)  -  ( 1  / 2tt i )  V1  V(z)(z-y)_1dz 
B'-i°° 


(21) 
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where  -1/2  <  0  <  Re  u  <  S'  <  0.  Inserting  the  ascending  series  of 
[33]  into  the  integrand  of  eqn.  (17)  and  integrating  term  by  term  re¬ 
sults  in 

a  _ 

/  K  (r»;  s/D)(l//r)dr  =  (it  t/a/2sinpn)  Z  (1/n! ) - 
0  y  n=0 

U2n"y/r(n-p+l)(2n-p+l/2)  -  £2n+u/r(n+y+l )(2n+y+l/2)}  (22) 

which  is  convergent  for  | Re  p |  <  1/2.  On  extending  the  range  of  y  to 
all  the  complex  values  by  analytic  continuation,  substituting  eqns.  (22) 
and  (17)  into  eqn.  (21)  and  shifting  the  path  of  integration  to  the 
left,  simple  poles  occur  at  z  =  -2n  -1/2  (n  =  0,  1 ,  2,  . . . )  with  resi¬ 
dues  (-l)n  £2n/  { (2n  +  u  +  1/2)  r(2n  +  1/2)}.  There  are  no  poles  at  z  * 
j,  (j  =  1,2,...)  because  of  the  relation 

oo  oo  _  _ 

Z  Z  (1/n!)  {rn+J/r(n+j+l)(2n+j+l/2)  -c  Zn'J/r(n-j+l)(2n-j+l/2)} 
n=0  j=l 

=  0  (23) 

Hence  Vp  becomes 

oo 

Vp  =  -(C°  *S/2s)  Z  (-l)n  C2n/r(2n+l/2)(2n+u+l/2)  (24) 

K  n=0 

The  function  W  can  be  expressed  as  a  sum  of  Wp  and  in  a  similar 
way.  If  eqn.  (19)  is  substituted  into  eqn.  (21)  with  W  and  Wp  replacing 
V  and  Vp,  respectively,  and  the  path  of  integration  is  shifted  to  the 
left,  simple  poles  arise  from  T(l/4  +  z/2)  at  z  -  -2n-l/2  (n  =  0,1,2,...) 
with  residues  (- 1 )n/n !  and  from  r(z)  at  z  =  -j  (j  B  1,2,...)  with  resi¬ 
dues  (- 1 ) J/ j I .  Consequently 
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WD  »  -  z  r(n+l/2)  4(2n+l/2)c4n+1/4n!r{2n+l/2)r(2n+3/2)(2n+p+l/2) 

K  n=0 

+  I  r(j/2+l/4)<Kj)C2j/4(j-l)!j!r(j/2+3/4)(j+y)  (25) 

3=1 

Although  the  integral  part  of  eqn.  (18)  is  analytic  for  all  values 

of  y,  X  is  analytic  in  the  restricted  region  Re  y  <  1/2,  owing  to  Lp(y)/ 

r(-y). 

Combining  eqn.  (20)  and  W  =  Wp  +  with  eqn.  (16)  yields 
VN  +  X  ‘WN  =  'VP  +  WP  +  LPi/2  =  G(v)) 

where  we  have  represented  both  sides  as  G(u).  Since  the  left  hand  side 
and  the  middle  are  analytic  at  least  in  the  regions  Re  y  <  0  and  Re  y 
>  -1/2,  G(y )  is  analytic  in  the  strip  -1/2  <  Re  y  <0.  However  we 
can  let  G  be  analytic  over  the  entire  region  of  y  by  analytic  continua¬ 
tion.  Since  £=  0(l/y),  Lp  =  0(1/41),  Vp  =  0(l/y)  and  Wp  =  0(l/y)  as  y 

approaches  infinity,  G  becomes  zero  as  y  approaches  infinity.  If 
Liouville's  theorem  [31]  is  applied  to  two  kinds  of  features  for  G,  G  is 
found  to  be  identically  zero.  Hence 

*(y)  =  2(Vp-Wp)/Lp  (26) 

Inserting  eqns.  (24)  and  (25)  into  eqn.  (26),  we  can  readily  express  4  as 

co 

*(u)  =  -{CVa-/4FsLp)  I  an(uKn  (27) 

n=0 

The  functions  an  are  functions  of  y  and  can  be  successively  determined 
to  give 


a0(u)  -  2/(2y+l ) 


(u)  =  4/ tt  (2  u+1 ) 

a2(y)  =  -8/3(2u+5)  +  8/tt2(2u+1  )  -2/3(v+l) 
a3(u)  =  16(1  /ir2  -  2/9 )/tt (2p+l )  -  4/3tt(p+1) 


The  total  current,  noting  eqn.  (4),  is  given  by 
_  a 

[TTnH  ■  -  /  2ttD(  3Cn(r,ir/2)/90)dr 
0  u 

=  -((2ttD)/(1+0)  /  v'(r,Ti/2)(l//r)dr  (28) 

0 


It  is  shown  in  Appendix  1  that  the  inverse  transform  of  eqn.  (13)  is 
given  by 

v'(r,TT/2)  =  (1/Tri )  /  (*(p)Cu/r(y))K  (r/s7ff)dy  (29) 

-ioo  p 


If  we  substitute  eqn.  (29)  into  eqn.  (28),  exchange  the  order  of  in¬ 
tegration  and  carry  out  integration  with  respect  to  r,  noting  eqn.  (22), 
we  have 


_  p  i°° 

[I/nF]  =  -(tt  D^a/in(l+s))  /  {<S(y)/r(y)sinyiO* 

-i« 

00 

I  (l/n!)U2n/r(n-p+l)(2n-u+l/2)  -C2n+w/r(n+y+l )(2n+y+l/2))dy 
n=0 


Completing  the  path  of  integration  in  the  right  half  plane  and  using 
eqn.  (23),  simple  poles  appear  at  y  *  2n+l/2.  When  the  residues  are 
calculated  and  eqn.  (27)  is  substituted  into  the  resulting  equation,  it 
follows  that 
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_  n  »  [k/2]  n  . 

[I/nF]  =  (4/b  CuDa/s(l+;))  I  i  (-1)  C  ak  ?n(2n+l/2)/T(2n+l/2) 

k=0  n=0  K'^n 

where  [k/2]  means  the  integer  part  of  k/2.  The  inverse  Laplace  trans¬ 
form  is  given  by 

I  =  (4nFC°Da/(l+;))[l  +  2/tt3/2/F  + 

Of  |( 

0/^)  :  I  (-l)k+n  r(k+l/2)  a?.  ?n+1(2n+l/2)T"k_1/2/r(2n+l/2)] 
k=l  n=0  CK-^n+i 

*(4nFC°Da/(l+c))  [1  +  0.35917  t-1/2  +  0.24648  t‘3/2 

+  0.20648  t‘5/2  +•••]  (30) 

? 

where  t  =  4Dt/a  . 

Small  Values  of  t  (Asymptotic  Expansion).  The  diffusion  equations 
in  circular  cylindrical  coordinates  are  given  by 

acyat  =  D  {(l/r)3(r3Cj/3r)/3r  +  32Cj/9z2}  (j  =  0  or  R)(31) 

The  initial  and  boundary  conditions  are  the  same  as  eqns.  (2),  (3),  (4) 
and  (5)  if  e  is  altered  to  z  and  z  is  set  to  be  zero.  On  making  a 
change  of  variable,  v  =  /r  (Cq  -  -£CR),  and  performing  Laplace  trans¬ 

formation,  eqn.  (31)  becomes 

sv/D  =  (l//r)3[r3(v//r)/3r]/3r  +  32v/3z2 

Applying  the  Hanckel  transformation  of  order  zero  [31],  given  by 

00 

U(p,z)  =  /  v(r,z)  A rp  JQ(rp)dr  (32) 


to  the  above  equation  yields 
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d2U/dz2  =  (p2  +  s/D  )U 

where  the  relation,  v  =  0(/r)  as  r  -*■  0,  has  been  used.  The  solution  of 

>2 

the  differential  equation  is  obviously  given  by  A*exp[-/p  -s/D  z]. 

When  we  eliminate  a  constant,  A,  from  U(p,0)  and  dll(p,0)/dz,  we  have 

/p2  +  s/D  U(p ,0)  +  dU(p,0)/dz  =  0 

If  we  rewrite  this  in  the  integral  form  with  the  aid  of  eqn.  (32)  as 
well  as  the  boundary  conditions  and  make  use  of  the  following  expression 
[35] 

a 

I  r  Jn(rp)dr  =  (a/p)  J,(ap)  (33) 

0 

then  we  obtain 

-C°a  J,(ap)/s/p  +  /  v(r,0)  /rp  JQ(rp)dr 
a 

a  _ 

+  M  /  v'(r,0)  /rp  Jn(rp)dr  =  0  (34) 

0  u 

where  M  =  M(p)  =  1//  +  s/D  and  v'(r,0)  =  3v  (r,0)/3z. 

We  introduce  three  new  functions: 

U+(p)  =  (1/2)  f  (v(r,0)  /rp  ^  (rp)/H, ^ ^ (ap)  /ap}dr  (35) 
a 

U_(p)  =  (1/2)  /{v(r,0)/rp  Hq^ (rp)/H] ^ (ap)/ap}dr  (36) 

d 

a  _ 

U°'(P)  =  /  v'(r,0)  /rp  JQ(rp)dr  (37) 

where  Hq^,  \  Hq^  and  H^2^  are  Bessel  functions  of  the  third 
kind  [36].  Inserting  the  asymptotic  expansions  of  Bessel  functions  of 
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the  third  kind  into  eqns.  (35)  and  (36)  and  letting  r  be  r+a  leads  to 

oc  . 

U.  (p)  %  +  (i/2)  /  7(r+a,0)  e^lprdr  (38) 

-  0 

for  large  values  of  Re  p. 

It  is  well  known  that  the  concentration  profile  can  be  represented 
by  the  profile  at  a  spherical  electrode  as  r  becomes  large.  Thus, 
v(r,0)  has  the  order  of  r'3,/2exp[-v/57E  (r-a)]  [25].  If  we  substitute 
this  into  eqn.  (38),  we  find  U+(p)  is  analytic  in  Im  p  >  -  Re  ATE 
and  U_(p)  is  analytic  in  Im  p  <  Re  /s/D  according  to  the  theory  of 
Laplace  transformation.  Expressing  eqn.  (34)  by  eqns.  (35),  (36)  and 
(37)  through  use  of  the  relation  2Jy  =  Hv^  +  H  dividing  the 
resulting  equation  by  H^^(ap)  v’ap  e71^4  (p  +  iATE)'1/2,  and  rear¬ 
ranging  it,  we  have 

U+/M+  +  UqM_/H-|  ^  }(ap)  Ap  +  R  +  S  =  0  (39) 

where 

M+  =  H+(p)  =  e"1'4  (p  +  i  A7E)'1/2 
M_  =  M_(p)  =  e_7Ti/4  (p-i  A7Ef 1/2 

R  =  -C°*^a/2spf1+  (40) 

S  MU.  -  C°A/2sp)Hj  ^(ap)/M+Hj  ^  ^ (ap)  (41) 

M+  and  M_  are  analytic  in  the  regions  Im  p  >  -Re  A/D  and  Im  p  <  Re 
ATE,  respectively,  and  their  product  equals  M. 

For  applying  the  Wiener-Hopf  technique,  R  is  expressed,  according 
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to  Cauchy's  integral  theorem  [31]  as 

R  =  R+  +  R_  (42) 

w+gi 

R+  =  (l/2ri )  f_K  +gi  (R(z)/(z-p))dz 

R.  =  -  ( 1 / 2r i )  /  (R(z)/(z-p))dz  (43) 

where  -Re  AJb  <  3  <  Im  p<  £'  <  0.  By  closing  the  path  of  the  integral 
in  eqn.  (43)  to  the  upper  half  plane  and  calculating  the  residue  at  z  * 
0,  R_  becomes 

R_  =  -(C°^a/2sp)(s/D)1/4 

Hence 

R+  »  R-  R.  =  (C°r'a/2sp)  {(s/D)1/4  -  1/M+}  (44) 

S  can  be  separated  into  two  parts. 


S  »  S+  +  S.  (45) 

S.  *  (l/27Ti)  761  (U  (z)  -  C°^a/2sz)  [H^tazJ/H/^aznUzMz-pJildz 
-°°+Bi  ‘  + 

(46) 

where  S_  corresponds  to  eqn.  (43).  Since  H^^(z)/H-j ^ \z)  %  -ie"2iz, 
the  singular  point  of  the  integrand  in  eqn.  (46)  is  restricted  only  to 
z  ■  -i  /s/D  coming  from  M+  in  the  lower  half  plane.  The  singularities 
in  the  bracket  of  eqn.  (46)  are  located  at  z  =  i  A 7U  and  z  *  0,  both 
of  which  are  pretty  far  from  z  =  -i  /s7U  for  large  values  of  |s|. 
Therefore  it  is  possible  to  let  z  be  -i  only  in  the  bracket  of 


eqn.  (46)  when  the  path  of  integration  is  shifted  to  the  lower  half 
plane.  Changing  the  variable  z  =  -i  /s 71)  +  iw,  carrying  out  integration 
and  employing  the  relation, 

U_(e'7rip)  =  -U+(p)  (47) 

we  have  the  following  leading  term: 

S+MU+(iA7&)  -  C°»W2is3/2}e'2a^7n/25/2  ^a3/2(iVs7U  +  p)  (48) 

On  combining  eqn.  (39)  with  eqns.  (42)  and  (45), 

R+  +  S+  +  U+/M+  =  -Uq  M./H-j  ^  ^(ap)v’ap  -  R_  -  S.  (49) 

R+  and  S+  as  well  as  U+/M+  are  analytic  in  the  region  of  Im  p  >  -Re 
v's/D,  while  each  term  on  the  right  hand  side  is  analytic  in  Im  P  <  0 
because 

Uq/h/^Up)^  'V  (1/2)  /  (iei(r_a)p  -  e-i(r+a)p)  v'(r,0)dr 

0 

Since  all  the  terms  on  the  left  hand  side  of  eqn.  (49)  tend  to  zero  as  p 
i°°,  it  follows,  according  to  Liouville's  theorem  [31], 

U+(P)  =  -M+(p)  (R+(p)  +  S+(p)}  (50) 

Substituting  eqns.  (44)  and  (48)  into  eqn.  (50)  and  letting  p  be  i*^7D, 
we  can  see  that  U+(i/s7U)  has  the  order  of  l/(sp).  Therefore,  S+  has 
the  order  of  (1/s)  exp  [- ✓T] ,  which  can  be  negligible  for  large  values 
of  | s | .  Then,  eqn.  (50)  becomes  asymptotically 

U+(p)  'v  (C°^a/2sp)  (l-M+(p)(s/D)1/4}  (51) 
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From  the  combination  of  eqn.  (47)  and  (51),  U_  is  given  by 
U_(p)  -v  (C°i/a/2sp)  {1  -  M_(p)(s/D)1/4} 

If  U+  and  l)_  are  eliminated  from  eqn.  (39)  with  eqns.  (40)  and  (41), 

U'q  takes  the  form 

Uq(p)  -v  (C°a/2Sv^p)(s/D)1/4  {h/^UpVM,  +  H1  (2)(ap)/M+}  (52) 

The  flux  for  the  total  current  is  given,  in  the  Laplace  transform, 
by 

_  a 

[1/nF]  =  /  2irrD  (9Cn(r,0)/3z)dr 
0  U 

a 

*  (2ttD/(1+;))  /  /r  v'(r,0)dr 
0 

Substituting  the  Hanckel  transform  of  U'(p,0),  given  by 

CO 

v'(r,0)  ■  /  lT(p,0)  v'rp  Jn(rp)dp 
0  u 

into  the  above  equation,  changing  the  order  of  integration  and  integrat¬ 
ing  with  respect  to  r  by  use  of  eqn.  (33),  we  have 

Co 

[I/nF]  =  (2TrDa/(l+c))  /  (J,  (ap)lT(p,0)/v'p)dp 

0  1 

U'(p,0)  is  equal  to  U'q(p)  since  v'(r,0)=  0  for  r  >  a  due  to  the  boundary 
condition  on  the  insulated  wall. 

Hence,  replacing  IT  by  U'q  (given  by  eqn.  (52))  leads  to 

[W]  *  (irC°Da2T/(l+Os)(s/D)1/4 


where 


T  =  /  (H] ( 1  ^(ap)/M_  +  H-j (^^(ap)/M+}  (J1  (ap)/p)dp  (53) 

The  integral  in  eqn.  (53)  is  evaluated  in  Appendix  2  to  give 

T  =  (s/D)1/4{1+(1//^)  I  r  (2n-3/2)r(n+l/2)2C1-2n/r(3/2-n)r(5/2-n)r(2n)2} 

n=l 

The  inverse  transform  for  the  total  current  is 
I  'o  (4nFC°0a/(l  +  ;))  {^F/2/r 

oo  _ 

+  (/5F/2)  I  r(2n-3/2)r(n+l/2)ZTn_1/r(n)r(2n)Zr(3/2-n)r(5/2-n)} 

n=l 

=  (4nFC°Da/(l  + :) ){Jr,/2\rr  +  tt/4  -  3ttt/210  -  315ttt2/221 - }  (54) 

DISCUSSION 

The  first  term  in  eqn.  (30)  represents  the  solution  for  the  steady 
state,  derived  from  potential  theory  [37,38],  while  the  leading  term  in 
eqn.  (54)  corresponds  to  the  well-known  Cottrell  equation.  They  express 
satisfactorily  the  two  kinds  of  limiting  behaviors  of  i-t  curves  which 
are  intuitively  expected. 

The  function,  f,  given  by 

I  =  {4nFC°Da/(l+;)}  f(i)  (55) 

is  calculated  numerically  from  eqns.  (30)  and  (54),  and  the  values  are 
plotted  against  x  in  Fig.  2.  The  curve  computed  from  eqn.  (30)  meets 
the  curve  evaluated  from  eqn.  (54)  in  the  domain  1.4  <  x  <  3.2.  This 
fact  indicates  that  eqns.  (30)  and  (54)  suffice  to  describe  chronoampero 
metric  curves  of  real  systems  with  good  accuracy. 

The  contribution  of  the  edge  effect  is  revealed  in  the  second  and 
the  subsequent  terms  of  eqn.  (54).  The  ratios  of  the  edge  contribution 
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terms  to  currents  due  to  the  Cottrell  term  are  less  than  5%,  1056  and  30% 

for  the  values  of  Dt/a^  less  than  8  x  10"^,  3.2  x  10*  ^  and  2.9  x  10*^, 

2 

respectively.  On  the  other  hand,  if  values  of  Dt/a  are  larger  than 
53,  3.5  and  0.65,  the  discrepancies  from  the  steady  state  are  within  5%, 

10%  and  30%  of  the  steady  state  current,  respectively. 

Some  studies  on  diffusion-controlled  currents  at  a  stationary  disk 
electrode  have  been  directed  toward  evaluation  of  the  term  succeeding 
the  Cottrell  term.  According  to  Soos  and  Lingane  [6],  the  current  ex¬ 
pression  can  be  expanded  in  the  form 

1  s  -rta^nFc'V  D/'t  {1  +  a i  J  Dt/ a  +  •••}  (56) 

where  the  coefficient  a^  is  a  constant  to  be  determined.  Some  values  of 
a^  reported  so  far  are  listed  in  Table  1  together  with  the  value  calcu¬ 
lated  from  eqn.  (54).  The  latter  value  is  very  close  to  the  simulated 
and  experimental  values  which  were  obtained  in  many  iterated  runs  by 
Kakihana  et  al.  [9,10]  and  experimental  values  by  Ito  et  al .  [8]. 

Dimensionless  i-t  curves  obtained  analytically  by  Soos  and  lingane 
[6]  and  computed  with  digital  simulation  by  Kakihana  et  al.  [9,10]  are 
shown  in  Fig.  2.  Since  the  expression  by  Soos  and  Lingane  is  just  the 
sum  of  the  Cottrell  term  and  the  steady  state  term,  it  does  not  ob¬ 
viously  describe  the  real  i-t  curve.  The  simulated  curve,  C,  is  over¬ 
lapped  on  curve  A  for  small  values  of  t  while  it  deviates  from  curve  A 
as  time  elapses.  This  is  a  natural  result  if  it  is  noted  that  digital 
simulation  frequently  involves  errors  for  a  number  of  Iterations  of 
computation  corresponding  to  long  electrolysis  time.  Most  experimental 

1 

i 

A 


chronoanperometric  curves  have  been  measured  in  the  range  t  <  0.3  [9], 
in  which  they  are  in  agreement  with  both  simulated  and  our  analytical 
curves  within  3%  error. 

The  clever  analytical  solution  of  Oldham  [14]  is  stated  to  be 
accurate  for  t  <  0.16.  In  fact  it  overestimates  the  value  of  f(r)  by 
only  0.7%  for  t  =  0.16.  The  approximation  remains  reasonably  accurate 
for  larger  values  of  t.  For  example,  for  t  =  1 ,  where  the  solution  of 
Soos  and  Lingane  overestimates  f (x )  by  24%,  the  corresponding  error  in 
Oldham's  solution  is  only  2.8%. 

Sarangapani  and  DeLevie  derived  an  analytical  expression  for  the  ac 
response  at  a  finite  disk  electrode  [11].  However,  eqn.  (15)  in  their 
paper  does  not  satisfy  the  boundary  condition  at  r  +  «.  If  we  carry  out 
the  integration  of  their  eqn.  (19)  in  a  fashion  similar  to  that  given  in 
Appendix  2  and  take  the  frequency  to  be  zero,  the  result  should  be 
identical  with  eqn.  (30)  or  (54).  However,  the  two  leading  terms  pro¬ 
vide  I  'v  -2.67  nFC^Da  for  the  steady  current  and  I  'v  -nFC°a2/D  (irt)"^2 
for  small  values  of  t  as  a  result  of  integrating  and  performing  inverse 
Laplace  transformation  of  their  expression.  Neither  of  these  results 
are  reasonable. 
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Table  1.  Magnitude  of  the  edge  effect  given  by  the  coefficient  of  the 

first  correction  term  to  the  Cottrell  equation.  The  coefficient, 
a-|,  is  defined  by  eqn.  (56). 


Authors 

Values  of  a-j 

Methods 

References 

Lingane 

2.12  +  0.11 

experimental 

5 

Soos  and 

Lingane 

r-  *  2.257.... 

analytical 

6 

Flanagan  and 

1.92  (published) 

7 

Marcoux 

1.79  (corrected) 

simulated 

9,10 

Kakihana  et  al. 

1.74  %  2.14 

experimental 

1.77  n,  1.98 

simulated 

9 

Dayton  et  al . 

2.16  +  0.35 
(at  carbon  paste) 

experimental 

12 

3.21  +  0.27 
(at  carbon  fiber) 

Ito,  Asakura 
and  Nobe 

1.77 

experimental 

8 

Heinze 

1.80-2.20 

simulated 

15 

This  work 


n  =  1.772... 


analytical 


APPENDIX  1 


This  appendix  shows  that  the  inverse  form  of  the  transformation  de 
fined  by  eqn.  (13)  is  given  by  eqn.  (29).  Let  us  prove  that  the  right 
hand  side  of  eqn.  (29)  in  which  4>  is  replaced  by  eqn.  (13),  i.e., 

(l/ni )  /  v  /  (v'l  (qt)/t)dt  K  (qr)dp  =  J  (A-l) 

-i«  0  y  u 

is  equal  to  v',  where  q  =  /s/D.  Substituting  the  relation  [33] 

K  (z)  =  (72sinpi)  {I_y(z)  -  I;j(z)}  (A-2) 


into  eqn.  (A-l)  yields 


J  *  (l/~i)  /  /  (Ttuv72tsinuTT)  {I  (qt)I  (qr)-I  (qt)I  (qr)}« 

-ioc  0  U  u  VJ  u 

dtdp  {A'3) 


If  u  is  replaced  by  -vi  in  the  first  term  of  the  bracket  and  eqn.  (A-2) 
is  employed  for  I_^(qt)  -  I^(qt),  eqn.  (A-3)  becomes 

J  =  (l/TTi )  /  uT(y)  I  (qr)dp  (A-4) 

where 

T(p)  =  /  (v'(r,  tt/2)  K  (qr)/r)dr  (A-5) 

0  y 

Since  V'  vanishes  for  r  >  a,  it  is  possible  to  extend  the  upper  limit  of 
the  integral  to  infinity.  Then  T  becomes  the  Kontorovich-Lebedev  trans- 
fcm  (see  Eqn.  8)  of  v'(r,  -n/2 ) .  Eqn.  (A-4 )  can  be  separated  into  the 
integrals  from  -i“  to  0  and  from  0  to  i°°.  Changing  variables  p  =  ix  in 
the  former  integral  and  p  =  -ix  in  the  latter  gives 


J  =  ( i /tt )  /  x  T(ix)  (Iix(qx)  -  I_ix(qx)}dx 

If  eqn.  (A-2)  is  inserted  into  the  above  equation,  J  is  given  by 

00 

J  =  (2/tt^)  /  x  sinh(Tix)  T(ix)K,  (qx)dx 

0  1X 

This  is  the  inverse  Kontorovich-Lebedev  transform  [30]  for  eqn.  (A-5). 
Therefore,  J  is  equal  to  7(r,  tt/2). 
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In  this  appendix,  the  integral  value  expressed  by  eqn.  (53)  is 
evaluated.  We  introduce  a  function, 

A  =  /{H^(ap)/N_(p)  +  H^(ap)/N+(p)}  (Jj(ap)7p)dp 

where  N+(p)  =  e 711/4  (p  +  i  /s7E)q,  N_(p)  =  e_7T1'/4  (p-iv/s7U)q  and  v  is  any 
real  constant.  It  is  evident  that  A  tends  to  eqn.  (53)  if  v  and  q 
approach  1  and  -1/2,  respectively.  Expressing  the  Bessel  functions  of 
the  third  kind  in  A  as  a  combination  of  Bessel  functions  of  the  first 
kind  yields  [36] 

cc 

A  =  (1/isinvTr)  /  {[J  (ap)  -  e'v"1  J  (ap)]/N  (p) 

0  v 

-[J_v(ap)  -  ew1  Jv(ap)]/N+(p)}  (J^apj/pjdp  (A-6) 

The  integral  representation  of  a  product  of  Bessel  functions  is  given  by 
[39] 

C+i°°  ?y+u+l 

J,(z)J  (z)  =  (l/27Ti )  /  (r(-x)r(2x+v+2)(z/2)"  ^‘/rfx+v+l) 

v  c-i® 

r(x+v+2)r(x+2)}dx  (A-7) 

where  -v/2  -  1  <  c  <  0. 

Replacing  the  products  of  the  Bessel  functions  in  eqn.  (A-6)  by 
eqn.  (A-7)  and  changing  the  order  of  integrals,  we  obtain 

A  =  - ( 1  /2tts i nv7T )  /  +1‘°{(f2-e"v1,if1 )  /"  (p2x_v/N  (p))dp 
c'-i»  0 

-(f2-e'm1f1)  /  (p2x'v/N+(p))dp}(a/2)2x‘v+1dx  (A-8) 

1  1  0 
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where 

f1  *  f'j (x ,v )  •  r(-x+v)  r(2x-v+2)/r(x+i)r(x+2)r(x-v+2) 
f2  *  f2(x,v)  =  r(-x)r(2x-v+2)/r(x-v+i)r(x-v+2)r(x+2) 

We  have  set  x  to  be  x-v  in  the  integral  representation  for  J-j Jv  In  which 
the  integrand  is  f -j .  Therefore  the  path  of  the  integral  for  f2  is 
located  in  the  left  half  plane  while  that  for  f-|  is  to  the  left  of  x  = 

1.  If  we  carry  out  the  integration  in  eqn.  (A-8)  with  respect  to  p,  we 
have  [40] 

/"  (p^x-v/N . (p) )dp  *  e±(*-W2-q/2*l/4)ji  (5/D)*-v/2-q/2+1/2. 

0  - 

r(2x-v+l)r(-2x+v+q-l)/r(q)  (A-9) 

Since  these  integrals  are  convergent  for  q  >  0  and  for  0  <  Re  x  <  Re 
(q/2) ,  q  is  temporarily  kept  positive  until  integration  is  carried  out 
with  respect  to  x.  After  that,  q  will  be  extended  to  all  real  values  by 
analytic  continuation  and  then  be  set  equal  to  -1/2. 

Substituting  eqn.  (A-9)  into  eqn.  (A-8)  and  rearranging  It  yields 

A  *  (2(s/D)'q/2/2TTir(q)sinv7r/  (a2s/4D)x_v/2+1/2. 

c'-i® 

T(2x-v+l  )r(-2x+v+q-l )  {f2sin(-x+v/2+q/2-l/4)it 
-fj  sin  (-x-v/2+q/2-l/4)ir}dx 

When  we  take  the  limit  as  v  -*•  1  by  use  of  L' Hospital's  theorem,  we  have 

n/9  9  c'+1® 

lim  A  -  -(2(s/D)'q/V2/ir(q))  / 
v+1  c'-i® 


(F2(x)+F1(x))dx 


(A-10) 
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F-j (x)  =  h.|  (x)f -j  (x ,  1 )  {h2(x)+o(l-x)-(Ti/2)cot(-x+q/2+l/4)Tt} 

F2(x)  «  h1  (x)f2(x»l )  {h2(x)+^(x)+(TT/2)cot(-x+q/2+l/4)Tr} 

h^x)  =  (a2s/4D)xr(2x)r(-2x+q)sin(-x+q/2+l/4h 

h2(x)  =  -(l/2)£n(a2s/4D)  -  ^(2x)  +  i^<q-2x)  -  i^(2x+l)  +  ^(x+1) 

where  v(x)  is  the  Psi  or  Digamma  function  [41].  As  described  before, 
the  path  of  integration  in  eqn.  ( A- 10 )  is  on  the  left  of  c'  =  1  for  Fi 
and  of  c'  *  0  for  F2«  Since  poles  do  not  arise  from  r(q-2x)  or 
i|/(q-2x)  under  the  conditions  of  eqn.  (A-9),  poles  resulting  from  closing 
the  contour  to  the  left  half  plane  are  the  simple  poles  at  x  *  0  coming 
from  three  terms  involved  in  the  parenthesis  of  F2  and  double  poles  at  x 
*  1/2  -n  (n  =  1 ,2,3, • • • )  coming  from  r(2x)^(2x+l ).  One  should  take 
notice  that  there  is  no  pole  at  x  =  -n  (n  *  1,2,3,***)* 

Each  contribution  of  simple  poles  at  x  =  0  can  be  calculated  in  the 
familiar  method  to  give  - { 1 /2 ) ( s/D) 1  ,  (s/D)1^4  and  (l/2)(s/D)^4  in 

turn,  after  taking  the  limit  q  -  -1/2.  If  we  let  the  residue  of  eqn. 

( A- 1 0 )  at  x  *  0  be  Aq,  we  obtain 

A0  «  ( s/D ) 1 /4 

In  order  to  calculate  the  residue  at  x  =  1/2-n,  which  is  symbolized  by 
An,  we  should  evaluate  a[{F2(x)+F1 (x)}(2x+2n-l )2]/3x  at  x  *  1/2-n.  If 
the  following  relations: 

4>(l/2-n)  *  ^(1/2+n) 

and 


x«1/2-n  " 


[3(iMx)  -  ipO-x)  +  itcot(-x+q/2+l/4)u)/3x] 


2tt2  (q  -  -1/2) 


are  used  for  the  calculation  of  the  residue,  An  can  readily  be  obtained 
An  *  (l/v€)(s/D)1/4(a2s/4D)1/2'nr(2n-3/2)r(n+l/2)2/r(3/2-n)» 

r(5/2-n)r(2n)2 

The  sum  of  An  provides 

oo 

T  *  lim  A  *  (s/D)1/4{l+(l//iF)  I  (a2s/4D)1/2‘n  r(2n-3/2)- 

v  -  1  nml 

q  "  "1/2  2  ? 

r(n+l/2)Vr(3/2-n)r(5/2-n)r(2n)2} 
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FIGURE  CAPTIONS 

Fig.  1.  Spherical  and  cylindrical  coordinates  located  on  the  stationary 
micro  disk  electrode.  Dashed  vector  for  the  cylindrical  coor¬ 
dinate;  solid  vector  for  the  spherical  coordinate. 

Fig.  2.  Dimensionless  chronoamperometric  curves. 

A:  calculated  from  eqn.  (54),  B:  calculated  from  eqn.  (30), 

C:  computed  by  means  of  digital  simulation  [9,10],  D:  curve  by 
Soos  and  Lingane  [6].  The  ordinate  is  the  function  f(t)  given 
by  eqn.  (55). 
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